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Abstract 

Recent data on the structure function F2{x, Q^) at small values of x are analysed 
and compared with theoretical expectations. It is shown that the observed rise 
at small x is consistent with a logarithmic increase, growing logarithmically also 
with . A stronger increase, which may be incompatible with unitarity when 
extrapolated to asymptotically small values of x, cannot be inferred from present 
data. 



1 Introduction 



The HI and ZEUS collaborations at the ep collider HERA have published measurements 
of the structure function F2{x, Q"^) which extend the range in x and by two 

orders of magnitude as compared to previous fixed target experiments. A prominent rise 
of the structure function has been observed at small values of x, which has stimulated a 
variety of theoretical investigations. 

The behaviour of the structure function F2{x,Q'^) for x —>■ 0, with fixed, cor- 
responds to the high-energy or Regge limit for virtual photon-proton scattering. The 
computation of the high-energy behaviour of cross sections in QCD is an important and 
still unsolved problem. Theoretical approaches to understand the Regge limit are almost 
exclusively based on perturbation theory. The standard evolution equations^ predict a 
rise of the structure function F2 at small x. This holds for the double-asymptotic solution 
of the evolution equations at large and small x [Q, as well as for the solution of the 
BFKL equation ^ which makes a prediction for the growth at small x for fixed Q^. One 
expects that the rise of cross sections at large energies is eventually damped by screening 
corrections leading to an asymptotic behaviour which, for proton-proton scattering, 
has to satisfy the Froissart bound 0. 

The starting point of this paper is the examination of the small-x regime of the 
HERA data. For values of x below about 10^^ all measurements are compatible with 
a double- logarithmic behaviour in x and Q^. We then address the question whether 
the observed rise of the structure function F2 at small x is consistent with unitarity 
bounds. It is argued that a In ^ increase, if persistent to asymptotically small x , may 
be compatible with constraints from unitarity. The double-logarithmic fit to the data is 
compared with the double-asymptotic form considered by Ball and Forte M, WO . 



2 Phenomenological analysis 

The collaborations HI and ZEUS have provided measurements at various values of 
covering the range in x as shown in fig. |l|. The measurements of both collaborations are 
compatible with each other. The 1994 data by the HI collaboration ^ are used for the 
fits, since at present they are the most precise ones, and they also allow the distinction 
between correlated and uncorrelated contributions to the systematic uncertainties. 
^For a review, see 0. 
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Figure 1: (Q'^ ,^)-phasespace region covered by the HI and ZEUS measurements of F2. 

The behaviour of the structure function F2 in the covered kinematical range is dis- 
played in figs. H and || as function of x and Q^, respectively. The Q^-dependence for 
three bands in x, displayed in fig. 0, shows a logarithmic behaviour. The slopes are 
strongly x -dependent. For x ^ 0.05 the slope vanishes and the structure function F2 
becomes independent of Q^. Fig. |^ shows the x -behaviour of F2 for three values of Q^. 
The striking new feature of the HERA data is the prominent rise at values of x below 
about 10~^. This small-x regime connects to the valence region, which was intensively 
investigated in previous low energy experiments. 

In the following only the small- a; data will be considered. The restricted phase space 
region is defined by 

X < 0.010 , >5 GeV^ . (1) 

The cut in x implies an effective cut of the high Q^-data (cf. fig. |l]). 

It is a remarkable property of the data that in this small-x regime they are well 
described by a polynomial linear in log ^ for every measured Q^-value, 

F2{x, Q') = uoiQ') + u^m log ^ . (2) 

X 

Possible higher order terms in log ^ are statistically not significant and are consequently 
not considered. The quantities Uq and Ui are functions of Q^, as well as the quantity logf , 
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Figure 2: F{x , Q"^) versus for three values of x : upper points 0.0001, middle points 
0.003, lower points 0.05. 



q' = 35 GeV^ 



= 4.8 GeV^ 



= 200 GeV^ 



Figure 3: F{x , Q^) versus x for three values of : 4.8, 35 and 200 GeV^. 
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which is uniquely defined as the weighted average (log ^). Using uncorrelated errors, also 
the uncertainties of uq and ui are uncorrelated. The numerical value of f (Q^) refiects the 
available x-range for any given (see fig. |l]), as well as the precision of the data. It is 
then possible to represent the whole body of data in the restricted phase space region, 
defined by (Q), for each measured Q^-value by three numbers, 

v{Q^) , uoiQ^) ± 6uo{Q^) , Mi(Q2) ± Sui{Q^) . (3) 

In terms of F2 the two independent functions uq and ui represent average and slope, 

= -^F2{x,Q'). (4) 
c^log- 

The data on F2 is consistent with a linear dependence in both logQ^ and log^. 
Furthermore, its extrapolation to smaller values in logQ^ and log ^, respectively, suggests 
the existence of a common "fixpoint" {xo,Qq) (cf. figs. |^, All this can be summarized 
in an ansatz for F2 which is linear in the double-logarithmic scaling variable ^, 



— log - 

Qo ^ 



^ = log 7^ log — • (5) 



This simple form, when confronted with the data given in (P), implies 



91] ^ 



a + mlog-^log— ^ = Mo(Q' 



'^logTil = «i(<3 ) ■ (6) 
The first of these equations can be cast into the form 

log [xq/v) 

thus allowing the comparison between the directly measured slope ui and the one ob- 
tained from the point (f , uq) and the fixpoint (xo,a). This is illustrated in fig. |^. The two 
parameter pairs (a,xo) and (m,QQ) are correlated, as is obvious from eq. (||). They are 
chosen as follows : 

a = 0.078 , m = 0.364 , xq = 0.074 , = 0.5 GeV^ . (8) 

Using the Hl-data with uncorrelated errors the measured structure function F2 is 
plotted in fig. ^ as a function of the scaling variable ^ computed from x and with the 
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Figure 4: The open circles represent the directly measured slopes ui, the full circles are 
the induced slopes for xq = 0.074 and a = 0.078. 




Figure 5: The F2-data of HI- 1994 is plotted versus the scaling variable^. The open circles 
represent the data for Q"^ above 5 GeV^; the crossed circles correspond to the data with 
Q'^ between 3 and 5 GeV^. 
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parameter values for xq and Ql as given in eq. (H). The x^/dof is 83/72. It turned out 
that even the data below 5 GeV^ are well described by eq. (||). We have checked that 
also the other data on F2 (cf. p, |l|) are well described by the same set of parameters. 



3 Constraints from unitarity 

No rigorous bound is known on the asymptotic behaviour of the structure function 
-^2(3^5 Q^) as X — 0, with fixed. However, within the framework of the parton model 
constraints on the allowed growth of parton densities at small x can be obtained by 
considering the Froissart bound on the total proton-proton cross section at large 
energies, 

atot{s)<^(\n-y . (9) 

Here 771^^ is the pion mass, s is the center-of-mass energy squared and Sq is an unknown 
constant. The total cross section measured at the Tevatron is of order it /mi ~ 100 mb. 

The cross section for the production of particles with high tranverse momentum, 
p5_ > A*^ ^ ^QCD5 can be calculated perturbatively. Clearly, this part of the cross 
section, ap^rt, must also satisfy the Froissart bound, i.e. 

Crpert{s,p']_ > /i^) < (7tot{s) . (10) 

Here apert is evaluated in terms of parton cross sections and parton densities. Eq. (p!0|) 
then yields a consistency condition for the behaviour of the parton densities at small x. 

In the parton model the dominant contribution to (Jpert is elastic gluon-gluon scat- 
tering. Consider the production of two gluons with transverse momentum p± ± Ap± and 
rapidities 1/3 ± Ay and Ay, respectively. Other partons in the final state are summed 
over. The corresponding cross section reads in lowest order perturbation theory (cf. 



AcTgg ~ Casifxfxigixi, n)x2g{x2, n)-^F{e*)Ap'iAy'^ . (11) 

P± 

Here Xi and X2 are the momentum fractions of the gluons in the initial state, and g{x, fi) is 
the gluon density at scale fi; C is a constant and F is the averaged squared matrix element 
of the gluon-gluon cross section, which depends only on 6*, the scattering angle in the 
gluon-gluon center-of-mass system. Note, that the multiplicity factor, which connects 
inclusive and exclusive cross sections, is (1 + 0{af)) in eq. (|TlD . In the special case 
Vs + 1/4. = 0, one has 

xi = X2 = X = . (12) 
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The cross section (|TTj) is a tiny fraction of the total cross section, Adgg < ~ l/xb. 

Its dependence on the center-of-mass energy s scales like {xg{x, fi))'^, with x = Ij)^! \fs- 
Clearly, if the gluon density satisfies the bound 

xg{x,fi) < B\n— , (13) 

X 

where 5 is a constant, then the corresponding bound on Aagg and the Froissart bound 
scale with s in the same way, i.e. like (Ins)^. Hence, Acr^g will always remain far below 
the Froissart bound. 

One may hope to derive a stronger bound on the gluon density based on a complete 
evaluation of the perturbative cross section cXpert- Already the integration over the ra- 
pidities of the two final state gluons yields another factor of (Ins)^. Contributions with 
additional gluons in the final state yield further powers of Ins which eventually build 
up the full BFKL ladder [^]. It is conceivable that these perturbative corrections can be 
absorbed in a properly defined gluon density for which the bound (0) may then apply. A 
similar analysis could be carried out for deep-inelastic scattering where one has only one 
gluon in the initial state. In this way it may be possible to obtain a true bound on the 
structure function F2 at small x. This, however, is beyond the scope of our paper. In any 
case, it is clear that with a logarithmic growth of the gluon density the perturbative cross 
section can be extrapolated several orders of magnitude beyond present center-of-mass 
energies before it possibly reaches the Froissart bound. 

One can readily evaluate the contribution of photon-gluon fusion to the structure 
function F2 for a gluon density saturating the bound (]T3|). One finds (cf. |I2|), 

AF^ix, Q')=A + ^j:elB\n%ln^, (14) 

OTT g C^Q X 

where the sum extends over all quarks with masses small compared to Q, and the con- 
stant A depends on the renormalization scheme. Eq. (|1^ is identical to the fit (|^) of F2, 
obtained by the phenomenological analysis in sect. 2, if parameters are properly iden- 
tified. /^From the slope m = 0.364 one obtains i? ~ 3. At small x the gluon density is 
large, and the structure function is likely to be dominated by photon-gluon fusion.^ 

Eq. (|T^) relates the behaviour of F2 at small x to a constraint on the gluon density 
obtained from the total proton-proton cross section. Our discussion has been based on 
a comparison of cross sections in perturbation theory to leading order, and the effect 
of higher order corrections is not clear. It is conceivable that, at fixed Q^, F2{x,Q'^) 
similar fit to F2, with In l/x replaced by a small power x^^ has been performed in ]l3|. 
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continues to rise like In ^ down to asymptotically small values of x. In fact, such a 
behaviour has been predicted by Bjorken. His starting point is an expression for the 



photon wave function renormalization constant derived by Gribov 

sR{s) 



a 
3^ 



ds- 



(15) 



+ s)2 ' 

where R{s) is the hadronic cross section in e^e~-annihilation normalised to the /x-pair 
cross section. Based on the aligned-jet picture of deep inelastic scattering he then obtains 



for the inclusive structure function 

F2{x,Q^) (X R \n 



1 



X 



(16) 



This corresponds to a gluon density saturating the bound {^^- At present, however, we 
do not know whether eqs. ( P^ ) or (|T^) represent the correct behaviour of F2{x,Q'^) at 
asymptotically small values of x. There is neither a proof that such a behaviour will ever 
be reached nor can it be excluded that it might set in already at moderate values of x. 



4 Comparison with double-asymptotic scaling 

So far we have discussed the Regge limit, i.e. a; with fixed. In this limit the double- 
logarithmic scaling form (^) of the structure function F2 may be correct. However, in an 
appropriate simultaneous limit Q"^ ^ 00 and x — (see below), the double- logarithmic 
scaling form is expected to be incorrect. In this limit, for appropriate boundary condi- 
tions, the structure function is given by an asymptotic solution of the standard evolution 
equations , which was found more than 20 years ago P| , and which has recently been 
thoroughly studied by Ball and Forte |P, |10|- This solution is known to correspond to a 
summation of all terms of the form 

(..ingmi)" . (17) 

Hence, at small x, it increases faster than any power of In ^. The double-logarithmic form 
(^) corresponds to the first term in this series. A comparison of these two expressions 
with data on F2 therefore tests for the evidence of terms more singular than In ^ at small 

X. 



The asymptotic solution is conveniently expressed in terms of the variables [|10 
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Figure 6: Comparison of F2 -data with the Ball-Forte fit (right vertical scale) and the 
double-logarithmic fit (left vertical scale); the horizontal axis corresponds to 2.2 a and 
^, respectively. 



where as is the two-loop QCD couphng, 



«^W) = 7r^ • (19) 



47r / A Inlnf ^ 

The first two coefficients of the /5-function and other relevant parameters read pOf 

2 S8 19 1 / 2 \ 

A = ll--n,, A = 102--„,. y = -. ^ = _(ll + -„,). 

1 / /5i 103 \ 78 



At large values of a the structure function F2 can now be written as |]T0| 

F2(a;,Q^) = C exp27(T exp f — 5— — - ln7(T — In — ] 

V P 2 7/ 

X (1 + (e+ + 6_) «,(Q) - 6+«,(go)) - f 1 + + . (21) 

7 V ^ P / 

Here C is an unknown normalization constant. For large values of a and p the ffist 
exponential gives the dominant contribution, 

F2(x,Q2) exp(27(T) . (22) 
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In this limit the structure function only depends on cr, which is referred to as double- 
asymptotic scaling. 

Eq. (^) gives F2 in terms of four a priori unknown constants, A, Qo? and C . 
It is known that a proper choice of these parameters yields a good description of the 
measured data on F2 in the range of small x and large [0, ^. Using the optimised 
values of the constancy of the ratio between the data and the theoretical prediction 
( PTD is clearly borne out in fig. ^ and quantified by "x^ jdof = 93/69. Eq. has been 
used for the number of flavours nj = A. A more sophisticated treatment incorporating 
Q^-dependent threshold effects has not been attempted. 
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Figure 7: The structure function F2 for = 15 GeV^ for the Ball-Forte fit (open circles) 
and the double-logarithmic fit (stars) extrapolated to values of x one decade below the 
HERA range. 

What does this mean concerning the validity of double-asymptotic scaling? As we 
have seen in sect. 2 (cf. also fig. |^), the data are equally well described by just the 
first term in the sum of double logarithms given in eq. (0). Hence, the characteristic 
feature of double-asymptotic scaling, a growth stronger than any power of In ^ cannot 
be inferred from the present HERA data. This more singular behaviour should become 
visible if, at given the range in x is extended by one order of magnitude. As fig. |^ 
illustrates, a clear difference is then expected between the double-logarithmic fit (H) and 
the double-asymptotic form 
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5 Conclusions 



We have shown that data on the structure function F2{x, Q"^) at small x and large 
published up to date are consistent with a logarithmic growth in ^ as well as Q^. This 
is of interest for several reasons. 

No rigorous bounds are known on the asymptotic behaviour of F2{x,Q'^) in the 
Regge limit, i.e. as x — > with fixed. However, according to our discussion in sect. 3, 
a logarithmic growth of F2 can be extrapolated to smaller values of x by many orders 
of magnitude without getting into conflict with the Proissart bound. Such a logarithmic 
increase with | may even be the correct asymptotic behaviour of the structure function. 

In the double-asymptotic regime of small x and large pcrturbative QCD predicts 
double-asymptotic scaling for the structure function F2, given sufficiently soft input 
distributions. This implies a growth stronger than any power of In^. So far, however, 
the data are still consistent with double-logarithmic scaling corresponding to an increase 
like In^. This may mean that at HERA the small-x regime has not yet been reached, 
where the strong growth expected on the basis of perturbative QCD will become visible. 

More precise measurements may differentiate between double-logarithmic scaling, 
double-asymptotic scaling and the even more singular BFKL power behaviour. As dis- 
cussed in sect. 4, these differences should become manifest if, for given Q^, the range in 
X is extended by one order of magnitude. This corresponds to an increase in the center- 
of-mass energy squared by one order of magnitude, which could be reached at future 
coUiders, such as LEP(8)LHC or at a 500 GeV Linear Collider (g) HERA. 

We would like to thank J. Bartels, P. V. Landshoff, P. Schrempp and P. M. Zerwas 
for valuable discussions. 
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